Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 2 February 2008 



(MN m&i style file v2.2) 



CMB polarization constraints on radiative feedback 

C. Burigana V L.A. Popa f, 2,1 R. Salvaterra |, 3 

R. Schneider §, 4 T. Roy Choudhury % 5 and A. Ferrara || 6 

1 INAF-IASF Bologna, Istituto di Astrofisica Spaziale e Fisica Cosmica di Bologna, 
Istituto Nazionale di Astrofisica, via Gobetti 101, 1-40129 Bologna - Italy 

2 Institute for Space Sciences, Bucharest- Magurele, Str. Atomostilor, 409, PoBox Mg-23, Ro-077125, Romania 

3 Dipartimento di Fisica G. Occhialini, Universitd degli Studi di Milano Bicocca, Piazza della Scienza 3, 1-20126 Milano, Italy 

4 INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 1-50125 Firenze, Italy 

5 Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK 

6 SISS A /International School for Advanced Studies, Via Beirut 4, 1-34100 Trieste, Italy 



2 February 2008 



ABSTRACT 

We compute the imprints left on the CMB by two cosmic reionization models 
consistent with current observations but characterized by alternative radia- 
tive feedback prescriptions {suppression and filtering) resulting in a different 
suppression of star formation in low-mass halos. The models imply different 
ionization and thermal histories and 21 cm background signals. The derived 
Comptonization, u, and free-free distortion, yB, parameters are below current 
observational limits for both models. However, the value of u ~ 1.69 x 10~ 7 
(~ 9.65 x 10 -8 ) for the suppression (filtering) model is in the detectability 
range of the next generation of CMB spectrum experiments. Through the 
dedicated Boltzmann code CMBFAST, modified to include the above ioniza- 
tion histories, we compute the CMB angular power spectrum (APS) of the TT, 
TE, and EE modes. For the EE mode the differences between these models 
are significantly larger than the cosmic and sampling variance over the multi- 
pole range I ~ 5 - 15, leaving a good chance of discriminating between these 
feedback mechanisms with forthcoming/future CMB polarization experiments. 
The main limitations come from foreground contamination: it should be sub- 
tracted at per cent level in terms of APS, a result potentially achievable by 
novel component separation techniques and mapping of Galactic foreground. 

Key words: Cosmology: cosmic microwave background - galaxies: formation 
- intergalactic medium. 



1 INTRODUCTION 

The accurate understanding of the ionization history 
of the universe plays a fundamental role in the mod- 
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crn cosmology. The classical theory of hydrogen re- 
combina t ion f o r pure baryonic cosmological m odels 
(|Peeblesl 1 19681 : IZel'dovich. Kurt fc Sunvaevl 1 19681 ) has 
been subseq uently extended to non-bar y onic dark mat- 
ter models |Zabotin Naselskvl Il982l : |j ones fc Wvsd 
ll985l : ISeager. Sasselov fc Scottll2000l 7 and recently accu- 
rately updated to include also helium recombination in 
the cu rrent cosmological scenario (see lSwitzer fc Hiratal 
(|2007l ') and references therein). Various models of the 
subsequent universe ionization history have been con- 
sidered to take into account additional sources of photon 
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and energy production, possibly associated to the early 
stages of structure and star formation, able to signifi- 
cantly increase the free electron fraction, x e , above the 
residual fraction (~ 10 -3 ) after the standard recombina- 
tion epoch at z ICC ~ 10 3 . These photon and energy pro- 
duction processes associated to this reionization phase 
may leave imprints in the cosmic microwave background 
(CMB) providing a crucial "integrated" information on 
the so-called dark and dawn ages, i.e. the epochs before 
or at the beginning the formation of first cosmological 
structures. For this reason, among the extraordinary re- 
sults achieved by the Wilkinson Microwave Anisotropy 
Probe (WMAP) mission, the contribution to the under- 
standing of the cosmological reionization process has 
received a great attention. 

To first approximation, the beginning of the reion- 
ization process is identified by the Thomson optical 
depth, r. The values of r compatible with WMAP 3yr 
data, possibly complemented with external data, are 
typically in the range ~ 0.06 — 0.12 (corresponding 
to a reionization redshift in the range ~ 8.5 — 13.5 
for a sudden reionization history), the exact interval 
depending on the consider ed cosmological mo del and 
combination of data sets l|Spergel et alj I2007T ) . While 
this simple "r-parametrization" of the reionization 
process and, in particular, of its imprints on the CMB 
anisotropy likely represents a sufficiently accurate 
modelling for the interpretation of current CMB data, 
a great attention has been recently posed on the 
accurate computation of the reionization signatures in 
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of WMAP accumulating data and of forthcoming and 
future experiments beyond WMAP. In this context, 
this w ork represents a step forward of our previous 
paper (|Schneider et al.ll2007l ) dedicated to the study of 
the impact of reionization, and the associated radiative 
feedback, on galaxy formation and of the corresponding 
detectable signatures. In that work we carried out a 
detailed comparison of two well defined alternative 
prescriptions (suppression and filtering) for the radia- 
tive feedback mechanism suppressing star formation 
in small-mass halos showing that they are consistent 
with a wide set of existing observational data but 
predict different 21 cm background signals accessible to 
future observations. We focus here specifically on the 
signatures detectable in the CMB. 

We observe that all viable (i.e. data-consistent) 
reionization models require radiative feedback. In fact, 
in the absence of radiative feedback Population (Pop) 
III stars would be allowed to form efficiently down to 
low redshifts, and, in order to match the low-z Gunn- 



Peterson data, a decreased Pop III star formation effi- 
ciency is required, in turn yielding a too low Thomson 
optical deptfQ. In other words, some feedback mecha- 
nism is necessary to have enough photons at high- 2 and 
to avoid at the same time an excess of photons at late 
times. Our paper indicates how to discriminate between 
the two most physical radiative feedback prescriptions 
with future CMB data. 

The paper is organized as follows. In Section [5] 
we briefly summarize the main theoretical and obser- 
vational aspects of the considered models focussing on 
those relevant for the analysis of the features in the 
CMB. In Section [3] we present the results of our com- 
putation for the CMB spectral distortions (Section l3.1[) 
and anisotropies in total intensity and polarization (Sec- 
tion I3.2|l and compare them with the foreground and 
sensitivity limitation of future CMB anisotropy missions 
(Section l3.3p . Finally, Section 2] summarizes our conclu- 
sions. 

Through this paper we assume a flat ACDM cosmo- 
logical model consistent with WMAP described by mat- 
ter and cosmological constant (or dark energy) density 
parameters Q m — 0.24 and J1a = 0.76, reduced Hubble 
constant h = f/o/(100km/s/Mpc) = 0.73, baryon den- 
sity Qbh 2 = 0.022, density contrast as = 0.74, and adia- 
batic scalar perturbations (without running) with spec- 
tral index n s = 0.95. We assume a CMB b ackground 
temperature of 2.725 K |Mather et al.lll999l) . 



2 REIONIZATION MODELS 

The latest analysis of Lya absorption in the spectra 
of the 19 highest redshift Sloan Digital Sky Survey 
(SDSS) quasars (QSOs) shows a strong evolution of the 
Gunn-Peterson Lya opacity at z ~ 6 (|Fan et al.ll2006l ; 
iGallerani. Choudhurv fc Ferrar a 2006). The downward 
revision of the electron scattering optical depth to r = 
0.09 ± 0.03 in the release of the 3-yr WMAP data 
(|Page et al.ll2007l ; lsp"ergel et al.ll2007l ). is consistent with 
"minimal reionization models" which do not require the 
presence of very massive (M > IOOMq) Pop III stars 
(|Choudhurv fc Ferraral [20061 ; iGnedin fc Fanll2006l ). We 
can then use these models to explore the effects of reion- 
ization on galaxy formation, to which we will refer to as 
"radiative feedback" . 

Recently, Choudhury & Ferrara (2005, 2006) de- 
veloped a semi-analytic model to jointly study cos- 
mic reionization and the thermal hist ory of the inter- 
;alact ic medium (IGM) . According to ISchneider et alj 
20071 ). the semi-analytical model developed by 



1 The filtering model implies a value of r lower than that 
obtained in the suppression model (as explained in the next 
Section); decreasing it further would push the value of r be- 
low the lcr limit set by WMAP 3yr data. 
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IChoudhurv fc Ferraral (|2005l ^ co mplemented by the ad- 
dition al physics introduced in IChoudhurv fc Ferraral 
l|200d ) involves: i) the treatment of IGM inhomo- 
geneities by adopting the procedure of Miralda-Escude, 
Haehnelt & Rees (2000); ii) the modelling of the IGM 
treated as a multiphase medium, following the thermal 
and ionization histories of neutral, HII, and Helll re- 
gions simultaneously in the presence of ionizing pho- 
ton sources represented by Pop III stars with a stan- 
dard Salpeter IMF ext ending in the range 1 — 100 Mq 
(|Schneider et alj|2006j ). Pop II stars with Z = O.2Z 
and Salpeter IMF, and QSOs (particularly relevant at 
2^6); Hi) the chemical feedback controlling the pro- 
longed transition from Pop III to Pop II sta rs in the 
merger-tree model by ISchneider et al.l l|2006l ); iv) as- 
sumptions on the escape fractions of ionizing photons, 
considered to be independent of the galaxy mass and 
redshift, but scaled to the amount of produced ioniz- 
ing photons. It then accounts for radiative feedback in- 
hibiting star formation in low-mass galaxies. This semi- 
analytical model is determined by only four free param- 
eters: the star formation efficiencies of Pop II and Pop 
III stars, a parameter, ?7 eS c, related to the escape frac- 
tion of ionizing photons emitted by Pop II and Pop III 
stars, and the normalization of the photon mean free 
path, Ao, set to repro duce low-redshift observatio ns of 
Lyman- limit systems (|Choudhurv fc Ferraral [20051 ) . 

A variety of feedback mechanisms can suppress star 
formation in mini-halos, i.e. halos with virial tempera- 
tures T V i r < 10 4 K ; particularly if their clusteri ng is 
taken into account ('K ramer. Haiman fe "Ohl|2006j). We 
therefore assume that stars can form in halos down 
to a virial temperature of 10 4 K, consistent with the 
interpretation of the 3-yr WMAP data (Haiman & 
Bryan 2006; but see also Alvarez et al. 2006). Even 
galaxies with virial temperature T v i r Si 10 4 K can be 
significantly affected by radiative feedback during the 
reionization process, as the increase in temperature of 
the cosmic gas can dramatically suppress their forma- 
tion. Based on cosmological simulations of reionization, 
IGnedinI (2000) developed an accurate characterization 
of the radiative feedback on low- mass galaxy. This study 
shows that the effect of photoionization is controlled by 
a single mass scale in both the linear and non-linear 
regimes. The gas fraction within dark matter halos at 
any given moment is fully specified by the current filter- 
ing mass, which directly corresponds to the length scale 
over which baryonic perturbations are smoothed in lin- 
ear theory. The results of this study provide a quanti- 
tative description of radiative feedback, independently 
of whether this is physically associated to photoevapo- 
rative flows or due to accretion supp ression. 

Following ISchneider et~ai1 (|2007l ). we consider here 
two specific alternative prescriptions for the radiative 

feedback by these halos: 

i) suppression model - following IChoudhurv fc Ferraral 



(|2006l l. we assume that in photoionized regions halos 
can form stars only if their circular velocity exceeds the 
critical value u C rit = 2fcflT/^m p ; here fi is the mean 
molecular weight, m p is the proton mass, and T is the 
average temperature of ionized regions, computed self- 
consistently from the multipha se IGM model; 
ii) filtering model - following IGnedinI (2000), we as- 
sume that the average baryonic mass Mb within ha- 
los in photoionized regions is a fraction of the uni- 
versal value fb — Qb/^im, given by the fitting for- 
mula Mb/M = f b /[l + (2 1/3 - 1)M C /M] 3 , where M 
is the total halo mass and Mq is the total mass of 
halos that on average retain 50% of their gas mass. 
A good approximation for the characteristic mass Mc 
is given by the linear-theory filtering mass, M^ 3 = 

(3/a) J a da'Mj /3 (a') l-(a'/a) 1/2 , where a is the 

cosmic scale factor, Mj = (47r/3)p (itci /Gpy'^ 2 , is the 
Jeans mass, p is the average total mass density of the 
Universe, and c s is the gas sound speed. 

The model free parameters are constrained by a 
wide range of observational data. ISchneider et all (|2007f ) 
reported the best fit choice of the above four parameters 
for these two models that well accomplish a wide set of 
astronomical observations, such as the redshift evolution 
of Lyman-limit absorption systems, the Gunn-Peterson 
and electron scattering optical depths, the cosmic star 
formation history, and number counts of high redshift 
sources in the NICMOS Hubble Ultre Deep Field. 

The two feedback prescriptions have a noticeable 
impact on the overall reionization history and the rel- 
ative contribution of different ionizing sources. In fact, 
although the two models predict similar global star for- 
mation histories dominated by Pop II stars, the Pop 
III star formation rates have markedly different red- 
shift evolution. Chemical feedback forces Pop III stars to 
live preferentially in the smalles t, quasi-unpolluted ha - 
los (virial temperature Si 10 4 K. ISchneider et al1 l2006'). 
which are those most affected by radiative feedback. In 
the suppression model, where star formation is totally 
suppressed below u cr it, Pop III stars disappear at z ~ 6; 
conversely, in the filtering model, where halos suffer a 
gradual reduction of the available gas mass, Pop III stars 
continue to form at z £ 6, though with a declining rate. 

Since the star formation and photoionization rate at 
these redshifts are observationally well constrained, the 
star formation efficiency and escape fraction of Pop III 
stars need to be lower in the filtering model in order to 
match the data. Therefore reionization starts at z £ 15 

in the filtering model and only 16% of the volume is 
reionized at z — 10 (while reionization starts at z ~ 20 
in the suppression model and it is 85% complete by 
z = 10). For 6 < z < 7, QSOs, Pop II and Pop III give 
a comparable contribution to the total photo-ionization 
rate in the filtering model, whereas in the suppression 
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Figure 1. Rcdshift dependence of the fraction of free elec- 
trons (top panel) and of the gas kinetic temperature (bottom 
panel) for the two considered models. See also the text. 



model reionization at z < 7 is driven primarily by QSOs, 
with a smaller contribution from Pop II stars only. 

The predicted free electron fraction and gas tem- 
perature evolution (see Fig. [T]) in the redshift range 
7 < z < 20 is very different for the two feedback mod- 
els. Bottom panel of Fig. [1] compares the evolution of 
the gas kinetic temperature and CMB temperature for 
the two models. In particular, in the filtering model the 
gas kinetic temperature is heated above the CMB value 
only at z 15. 

The Thomson optical depth, r = J Xen e aTcdt, can 
be directly computed for the assumed ACDM cosmo- 
logical model parameters given the ionization histories 
shown in the top panel of Fig.Q] We find r ~ 0.1017 and 
r ~ 0.0631 for the suppression and the filtering model, 
respectively. Note that these values are consistent with 
the Tho mson optical depth r ange derived from WMAP 
3yr data |Spergel et al.lr2~007l ) but with ~ la difference 
among the two models, leaving a chance of probing them 
with forthcoming CMB anisotropy experiments. 



3 SIGNATURES ON THE CMB 

The cosmological reionization leaves imprints on the 
CMB depending on the (coupled) ionization and ther- 
mal history. They can be divided in three categoriefl: i) 



2 Inhomogeneous reionization also produces CMB secondary 
anisotropics that dominate over the primary CMB compo- 
nent for I > 4000 and can be detected by upcoming exper- 

iments, like the Atacama Cosmology Telescope or ALMA 
l lSalvaterra et al.ll2005l; llliev et alj|2007h . 



generation of CMB Comptonization and free-free spec- 
tral distortions associated to the IGM electron temper- 
ature incresase during the reionization epoch, ii) sup- 
pression of CMB temperature anisotropes at large mul- 
tipoles, i, due to photon diffusion, and in) increasing 
of the power of CMB polarization and temperature- 
polarization cross-correlation anisotropy at various mul- 
tipole ranges, mainly depending on the reionization 
epoch, because of the delay of the effective last scat- 
tering surface. 

The imprints on CMB anisotropies are mainly de- 
pendent on the ionization history while CMB spectral 
distortions strongly depend also on the thermal history. 

It is interesting to compare the imprints left in the 
CMB for the two models considered here to understand 
if they could be distinguished by forthcoming and future 
experiments. 

3.1 Spectral distortions 

The Compton scattering of CMB photons with 
the electrons heated during the reionization pro- 
cess implies the gene ration of a global Comp- 



tonization distortion (Zel'dovich fc Sunvaevl 1 19691 ; 
Zcl 'dovich. Illarionov fc Sunvaevl 1 19721 ) characterized 
by the Comptonization parameter it ~ (l/4)Ae/ej, 
where Ae/e^ is the fractional amount of energy ex- 
changed between matter and radiation. In addition, 
we expect also the generation of a free-free distor- 
tion, physically coupled to the Comptonization one, 
because of the bremsstrahlun g photon production 
proce ss (jKarsaz fc Latterl Il96ll ; iRvbicki fc LightmarJ 
Il979l ) in the hot gas. Its amplitude is characterized 
by the so-called free- free distortion parameter, y B - 
For late processes, as in this case, the resulting dis- 
torted spectrum can be fully described to a very good 
precision with the analytical form alism described in 
iBurigana. De Zotti fc Danesei (|l995l ) based on the com- 
putation of the Comptonization a nd free-free distortion 
parameters, by simply modifying ijBurigana et al. . 2004) 
the cosmic expansion time to take into account the 
cosmological constant (or dark energy) contribution, 
dominant at low redshifts. By exploiting the ioniza- 
tion and thermal histories shown in Fig. [1] we find 
u~ 1.69 xl0~ 7 , y B ~9.01xl0" 10 and it ~ 9.65 x 10" 8 , 
ys — 5.24 x 10 -10 respectively for the suppression 
and the filtering model: these val ues are clearly well 
below the CQBE/FIRAS l imits (jFixsen et all 1 19961 : 
ISalvaterra fc Buriganal [2002). The two models show 
similar ionization and thermal histories at z £ 6 while 

important differences are predicted at z £ 6. These 
explain the different spectral distortion levels generated 
in the two cases. In fact, considering only redshifts 
z > 6, we find u ~ 7.98 x 10" 8 , y B ^ 4.37 x 10~ 10 

and it ~ 1.05 x 10~ 8 , y B ^ 8.26 x 10~ n respectively 
for the suppression and the filtering model. We note 
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that, while the free- free distortion levels predicted in 
these models are too small to be detected even by long 
wavelength (i.e. at A ^ 1 cm) spectrum experiments 
with precision comparable to - or even significantly 
better than - that o f COB E/FIRAS (jKogutJ Il99rj ; 
iBurigana fc Salvaterral I2003T I. such Comptonization 
distortion prediction levels are comparable to those that 
could be in principle observed by a future generation of 
CMB spectrum experiments, in particular at millimetre 
wavelengths, able to improv e by m 30 — 100 times 
the COBE/FIRAS results ilFixsen fc Mather! |2002| ; 
IBurigana. Salvaterra fc Zizzol 12004 ) . In this perspec- 
tive, future accurate m easurements at long wavelengths 
(|Burigana et al.l I2007T ) could significantly improve the 
reliability of the detection of such Comptonization 
distortions, removing the approximate degeneracy 
in t he joint determination of early (Bose- Einstein 
like, ISunvaev fc Zel'dovichl (|l970h ) and late distortion 
parameters that remains in the presence of accurate 
spectrum measures only at millimetre wavelengths. 

3.2 Anisotropies 

The detailed computation of the reionization imprint in 
the CMB anisotropy angular power spectrum (APS), 
Ci, requires the use of dedicated Boltzmann codes, 
properly implemented to have the possibility of intro- 
ducing the adopted ionization history. We have mod- 
ified the public version (4.5.1 ) of th e CMBFAST cod<0 
(see e. g. ISeliak fc Zaldaxriagal (jl996h ) to properly re- 
place (|Popa et al. 20051 ) the simple step function (or 
Heaviside function) approximation for the ionization 
fraction adopted there to model the reionization pro- 
cess with the considered ionization histories, shown 
in Fig. [1] Our results are reported in Fig. [2] Hav- 
ing neglected for simplicity tensor perturbations, we 
focus here on the TT (total intensity, i.e. tempera- 
ture), TE (temperature-polarization cross-correlation), 
and EE polarization modes of the CMB anisotropy APS. 
In Fig. [2] we display also the APS of the foreground in 
the V band (centred at 61 GHz) of WMAP 3yr dat£0, 
a frequency channel where the foreground is found to 
be minimum (or almost minimum) in both temperature 
and polarization l|Bennett et al.ll2003l : IPage et af1l2007l ) 
at the angular scales larger than ~ 1° of relevance in this 
context. It is simply deri ved with the anafa st facility 
of the HEALPix package ijCorski et al.ll2005l ) analysing 
the difference between the original map and the CMB 
map and adopting the mask used by the WMAP team 
for the polarization analysis (covering ~ 74 per cent of 
the sky). For sake of simplicity, no attempt is made in 



3 http://www.cmbfast.org/ 

4 http:/ /lambda. gsfc.nasa.gov/product/map/current/ 

5 http://healpix.jpl.nasa.gov/ 
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Figure 2. APS of CMB anisotropies for the three consid- 
ered modes TT, TE, EE, reported in each panel and com- 
puted by introducing the ionization and thermal histories 
shown in Fig. [T] in the CMBFAST code (solid lines: suppression 
model; dashes: filtering model). Thick lines denote correla- 
tion, while thin lines denote anticorrelation (appearing for 
the TE mode at I > 50). The APS of the foreground, dom- 
inated by the Galactic contribution, is reported for compar- 
ison (dot-dashes). The WMAP 3yr V band data are consid- 
ered here. See also the text. [Results expressed in terms of 
thermodynamic temperature fluctuations]. 

these foreground plots to subtract the noise contribu- 
tion to the overall AP£0. We are in fact interested here 
to low multipoles (£ ;$ few tens) where the differences 
between these two models show up, as expected for rel- 
atively late reionization processes, and the signal power 
dominates over the noise one. Note that the relative dif- 
ference between the APS computed for the suppression 
and filtering model as well as the foreground impact 
compared to the CMB signal increase going from the 
TT, to the TE, and then to the EE mode. 



6 The relevance of the noise term - important at high mul- 
tipoles — obvioulsy increases going from the TT, to the TE, 
and then to the EE mode because of the decreasing of the 
signal-to-noise ratio. 
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3.3 Foreground and perspectives from Planck 
and future polarization experiments 

To understand if these classes of models can be distin- 
guished with forthcoming and future CMB anisotropy 
experiments we compare their relative difference with 
the cosmic and sampling variance limitation, ACt/Ce, 
(|Knox!l 1995b . i.e. neglecting the contribution by the in- 
strumental noise, and a potential residual foreground 
contamination. At the low multipoles (I ^ 20 — 30) rel- 
evant here (see Fig. [5} the direct (ideal) instrumental 
noise limitation is in fact not critical for current and 
future space missions. On the contrary, the experiment 
sensitivity and frequency coverage and, in particular, 
the rejection (by instrument design and/or by subse- 
quent data analysis) of the various classes of system- 
atic effects are critical to have the possibility to accu- 
rately subtract the foreground contribution and to ef- 
fect h^ly_achjcve_the ultim ate c osmic varian c e lim it (see 
e.g. iMennella etafl (|2004 ) and iPopa et~ai1 (120071 ). and 
references therein, for reviews devoted to these aspects 
in the context of the Planck mission). So, the results re- 
ported here apply to the case of the forthcoming Planck 
mission^, at least in the best case of an effective excel- 
lent rejection of all the potential systematic effects, and 
to the next generation of CMB polarization anisotropy 
experiments with a (nearly) all-sky coverage studied in 
the NASA (Inflation Probe, CMBPol) and the ESA (B- 
PoQ) context. 

Following Nasel skv fc Chia ng (2004), we compute 
thequantity S M ' N = {Cf -Cf )/[0.5(C £ M +Cf )] where 
the index M or N specifies the APS calculated for a 
given model (for simplicity, we omit here the indices 
TT, TE or TE; they are reported in each figure panel). 
A potential foreground residual contamination, result- 
ing from a non-perfect (possibly biased) component sep- 
aration, is modelled here as & Fore = fCf ore /[0.5((7f + 
Cf)), where C[ ore is the foreground APS and the ef- 
fective parameter / characterizes the accuracy of the 
component separation method in the considered range 
of multipoles. 

Figs. [3] and [4] summarize our results. We display 
5 M ' N for M = S (suppression model) and N — F (fil- 
tering model). In the case of the EE mode, we report 
also the results obtained exploiting the Heaviside func- 
tion approximation of the two reionization histories with 
corresponding optical depths, M = SH and N = FH. 
It is also interesting to consider 5 ' for M — S and 
N = SH and for M = F and N = FH in order to 
quantify the relative error implied by the step function 
approximation of a more complex time-dependent reion- 
ization history. 

As evident from Fig. [3] (top panel), the differ- 



7 http://www.rssd.esa.int/planck 

8 http:/ /www. th.u-psud.fr/bpol/index.php 
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Figure 3. Relative difference, S S ' F , between the (TT and 
TE mode) APS of CMB anisotropics for the suppression and 
filtering models reported in Fig. [2] (thick solid lines) com- 
pared with the cosmic and sampling variance limitation cor- 
responding to a sky coverage of ~ 74 per cent (region between 
the dotted lines). We report for comparison the APS from a 
potential residual foreground (dot-dashes) corresponding to 
two different values of /: 0.1 (dashed line) and 0.01 (three 
dots-dashed line). See also the text. 

ence between the two models is overwhelmed by the 
cosmic and sampling variance in the case of temper- 
ature anisotropy measures. The difference in the TE 
mode (see bottom panel of Fig. [3} turns to be well 
above the foreground limitation even for a removal ac- 
curacy of about 10 per cent. Unfortunately, it is above 
the cosmic and sampling variance limitation by a fac- 
tor of ~ 1.5 only and for a small range of multipoles 
(e ~ 7 - io0 On the contrary, for the polarization 
EE mode APS the difference between the two consid- 
ered models is significantly larger than the cosmic and 
sampling variance over a interesting range of multipoles 
(£ ~ 5 — 15), as shown by Fig. [4] (top panel). In this 
case the main limitation derives from a possible resid- 
ual foreground contamination: as evident, a foreground 
removal accuracy at per cent level (or better) in terms 
of APS is necessary to accurately exploit the informa- 
tion contained in CMB polarization about the cosmo- 
logical reionization process. Note that, possibly except 
for some sky areas at very low foreground contamina- 

9 The apparently large features at t ~ 50 result from the 
very small absolute value of the TE mode close to its change 
of sign with the consequence of large relative differences be- 
tween the two models but only on a very small multipole 
range (see also Fig. [2}. 
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Figure 4. Relative difference, S S ' F , between the (EE mode) 
APS of CMB anisotropics for the suppression and filtering 
models reported in Fig. [2] (thick solid lines) compared with 
the cosmic and sampling variance limitation corresponding 
to a sky coverage of ~ 74 per cent (region between the dot- 
ted lines). In the top panel we report also for comparison 
the APS from a potential residual foreground (dot-dashes) 
corresponding to different values of /: 0.1 (dashed line), 0.03 
(dots-dashed line), 0.01 (three dots-dashed line), and 0.003 
(long dashes). In the bottom panel we show also S SH ' FH 
computed by replacing the two considered models with those 
obtained assuming the simple step function implemented in 
the public CMBFAST code for the reionization history (thin 
solid lines) for the corresponding values of optical depth. The 
relative difference between the results based on the suppres- 
sion (resp. filtering) model and its step function approxima- 
tion with r ~ 0.1017 (resp. with r ~ 0.0631) is shown by the 
dashed line (resp. three dots-dashed line). See also the text. 



tion level (jCarretti. Bernardi fc Cortiglionil l2006). sim- 
ilar (resp. or better) accuracies are needed for the 
observation of BB polarization mode of the CMB 
anisotrop-J^l for ratios of tensor to scalar perturbations 
^ some x 10 -3 on typical sky regions at low/moderate 

foreground contamina tion level (resp. at low multipoles) 
(|La Porta et al.ll2006l ). This calls for a further progress 
in com ponent separation techn i ques in pola r izatio n 
(see e.g. iBaccigahipi et~aH d2004h: lstivoli et ail ll200rJ ): 
IPaee et all l|2007l ): lAumont fc Macias-Perezl (|2007h for 
promising developments) and for an accurate mapping 
of the (mainly Galactic) polarized foregrounds in radio 
and far-IR band£3 to complement microwave surveys 



10 http://www.b-pol.org/pdf/BPOL_Proposal.pdf 

11 See e.g. the talks by J. P. Leahy and J. -P. Bernard at 



http: / /www.th.u-psud.fr/bpol/talks_orsay .meeting/Talks. html 



by improving both the foreground physical modelling 
and the component separation accuracy. 

Finally, we observe that, particularly at I ~ 10 — 20, 
the step function approximation of the reionization his- 
tory adopted in the simple r-parametrization, although 
adequate to reveal the bulk of the difference between 
these models, implies an error on the EE mode compa- 
rable to the cosmic and sampling variancj^l (see bot- 
tom panel of Fig. [3} . Therefore, a proper treatment of 
the cosmological reionization history in the Boltzmann 
codes is necessary for a fine comparison with the high 
accuracy data from forthcoming and future CMB po- 
larization anisotropy experiments. 



4 CONCLUSION 

We have analyzed the implications for the CMB of 
the self-consistent semi-analytical model developed by 
Choudhury & Ferrara (2005, 2006) to describe the effect 
of reionization and its associated radiative feedback on 
galaxy formation, modelling the suppression of star for- 
mation in low-mass galaxies due to the increase in tem- 
perature of the cosmic gas in ionized regions according 
to two different feedback prescriptions, the suppression 
and filtering model. These two models predict different 
ionization and thermal histories, in particular at z ^ 6, 
resulting into Thomson optical depths r ~ 0.1017 and 
r ~ 0.0631, respectively, consistent with the WMAP 
3yr data. 

According to the usual semi-analytical approach 
to model CMB spectral distortions, we computed the 
Comptonization and free-free distortion parameters, u 
and j/s- Their are found to be well below the current 
observational limits. In particular, the values obtained 
for the Comptonization parameter (u ~ 1.69 x 10~ 7 
and ~ 9.65 x 10 -8 for the suppression and the filtering 
model, respectively) are in range accessible to a future 
generation of CMB spectrum experiments, in particular 
at millimetre wavelengths, able to improve by « 30—100 
times the COBE/FIRAS results. 

We have modified the Boltzmann code CMBFAST 
to introduce the considered ionization histories and 
compute the CMB angular power spectrum (namely, 
the APS of temperature (TT), temperature-polarization 
correlation (TE), and polarization EE modes) improv- 
ing the step function approximation adopted in the sim- 
ple "r-parametrization". The differences between the 
predictions of the two models are negligible for the TT 
mode and small for the TE mode when compared to 
cosmic and sampling variance limitation but are of par- 
ticular interest for the EE polarization mode. They are 

12 A similar analysis for the TT (resp. TE) mode shows that 
this kind of error is overwhelmed by (resp. about one order 
of magnitude less than) the cosmic and sampling variance. 
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significantly larger than the cosmic and sampling vari- 
ance over the multipole range I ~ 5 — 15, leaving a 
good chance of discriminating between the two feedback 
mechanisms with forthcoming and future CMB polar- 
ization anisotropy experiments, and in particular in the 
view of the forthcoming ESA Planck satellite that will 
launched in about one year. The main limitation derives 
from foreground contamination: it should be subtracted 
at per cent level in terms of APS, a result achievable 
with joint efforts in component separation techniques 
and in the mapping and modelling of the Galactic fore- 
ground. Also, a particular care should be taken in the 
control of potential systematic effects down to negligi- 
ble levels and in the precise reconstruction of the APS 
from the data at low multipoles. 

Our analysis indicates that forthcoming CMB 
anisotrop y polarization data t ogether with future 21 
cm data (Schneider et al 1 120071 1 will play a crucial role 
in breaking current degeneracies and constraining the 
cosmological reionization process and the astrophysical 
properties of the sources responsible for it at high red- 
shifts (z 6). 
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